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The feasibility of tetraquark detection is studied. For the ccMd tetraquark we show that in present 
(SELEX, Tevatron, RHIC) and future facilities (LHCb, ALICE) the production rate is promising 
and we propose some detectable decay channels. 
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I. INTRODUCTION 

The purpose of this paper is to assess the possibil- 
ity of detecting certain tetraquarks in present and future 
facilities. Among many possible tetraquarks, the dou- 
ble charm tetraquark Tec = ccud = DD* with quantum 
numbers IS^ = 01+ is particularly interesting since it is 
very sensitive to the chosen effective interaction: 

• It is very delicate, it is either weakly bound or 
slightly unbound with respect to the two-body 
hadronic decay D + D* . 

• Its structure can be either predominantly "molecu- 
lar" or predominantly "atomic" with consequences 
for the production and decay. 

Double charm tetraquarks were intensively studied by 
many authors. Various approaches were applied, from 
lattice QCD and chiral heavy quark effective theories to 
nonrelativistic potential models. It was shown that, al- 
though the predictions of these theories agree in the bar- 
ion and meson sector, they give dramatically different re- 
sults for tetraquarks. For this reason, the double charm 
tetraquarks present an important laboratory for discrim- 
inating between different hadronic models. 

Moreover, our estimates for the production cross sec- 
tion of such states gives us some hope that they can be 
experimentally detected in near future. 

The most important ingredient in the production of 
the Tec tetraquark is double charm production. Experi- 
mental data for such events are very puzzling. The pro- 
duction of prompt J/ip at B factories as well as the pro- 
duction of Sec at SELEX are much larger than expected. 
Therefore the comparison of Tec production with J/ipcc 
at B factories and the comparison Tec production with 
Sec at SELEX could shed some new light on the mecha- 
nism responsible for such large double charm production. 

In Section 2 we presents our results for the Tec pro- 
duction at high energy colliders were we believe that the 
dominant mechanism for the initial double charm pro- 
duction is double gluon fusion. From experimental data 
we also estimate phenomenologically the production of 
the Tec tetraquark at B factories and at SELEX. The 
results of detailed four body calculations in nonrelativis- 
tic constituent quark model (Section 3) encourage us to 



further investigate this state. Since we found the Tee 
tetraquark to be weakly bound, we propose (Section 4) 
the branching ratio between hadronic and radiative de- 
cays as the most promising mechanism for the detection 
of these states. 



II. PRODUCTION OF DOUBLE CHARM AT 
VARIOUS FACILITIES 

The most promising mechanism for the production of 
the Tec tetraquark is the formation of the cc diquark fol- 
lowed by hadronisation into ccud. An alternative mecha- 
nism would exploit binding of D and D* mesons if when 
they are produced with small relative momenta. One 
might expect that the latter mechanism could drastically 
enlarge the production rate if the dominant configura- 
tion is molecular. Due to the very messy environment 
in hadron colliders, however, such a weakly bound sys- 
tem would too soon dissociate into free mesons by the 
interaction with surrounding partons of initial hadrons. 

The first step is to create two cc pairs with the c quarks 
close in the phase space and in colour antisymmetric 
state, so that in the second step they bind into the cc 
diquark. The binding energy of such a system is ~ 200 
MeV [1]. In the third step the diquark gets dressed either 
with a light u or d quark into a ecu or ccd baryon or with a 
light ud antidiquark into the ccud tetraquark. The prob- 
abilities for this two types of dressing can be estimated 
using the analogy of a single heavy quark fragmentation. 
The branching ratio of the b ^ B and b ^ Ab produc- 
tion at the Fermilab and at LEP experiment is 0.9 and 
0.1 [12], respectively, therefore we expect the same ratio 
in the hadronisation cc — > Sec and cc —^ Tec, respectively. 

The double charmed baryons were probably detected 
at SELEX [11]. It was estimated that 40% of the singly 
charmed baryons they see result from the decay of dou- 
bly charmed baryons. The most probable mechanism for 
the double charm production at SELEX is production of 
the single cc pair in the processes gg cc or qq — > cc 
while the second cc pair is created in the fragmentation 
of the heavy quark c — > ccc. However, theoretically it 
is still unclear why the SELEX has such a large cross 
section for double charm production. Since SELEX is 
a fixed target experiment the cc diquark is most likely 
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to be produced with high lab momenta which might be 
helpful in the detection as discussed in [2] . But since SE- 
LEX found, with their cuts, only about fifty candidates 
for double charmed baryons, the statistics for detecting 
double charmed tetraquarks should be improved. 

Next, we look at the production and detection of the 
Tec tetraquark in B-factories. Since the total mass of 
four D mesons is close to the cm. energy, the c quarks 
created in this process have small relative momenta which 
is very important in T^c production. This feature also 
ensures a smaller number of additional pions created in 
the e+e" annihilation and thus a cleaner reconstruction 
of Tec- Belle [13], [14] has reported a measurement of 
prompt J/ip production in e+e^ annihilation at ^/s = 
10.6 GeV and found that the most of the observed J/tjj 
production is due to the double cc production 
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which correspond to [13], [14] 

a(e+e- ^ J/V'cc) = CSTtg-jg ± 0.17 pb 

or to about 2000 events from their 46.2 fb~^ data sample. 
The theoretical nonrelativistic QCD prediction for this 
process is an order of magnitude smaller [3] , [4] , [5] , [6] , 
so this process is still not well understood [7]. But it 
is very likely that the analogous mechanism would also 
enlarge the cross section for the prompt production of cc 
diquark and thus the Sec baryon and T^c tetraquark. The 
cc diquark production cross section can then be estimated 
to be ~ 0.15 pb, which correspond to ~ 10^ Sec [8] and 
about ^ 10"* Tec per year. 

We now present our calculation for the Tec produc- 
tion at high energy colliders. The two colliding nucleons 
in TeV machines can be considered as two packages of 
virtual gluons whose number is huge for low Bjorken-x. 
Therefore we expect [9], [10] that in these facilities the 
dominant mechanism for double charm production would 
be a double gluon-gluon fusion: (g + g) + (ff + (?) 
(c + c) + (c + c). 

The usual hard production mechanism is heavy quark 
production followed by fragmentation, however this 
mechanism docs not include all the possible Feynman 
diagrams. In ref. [9], [10] it has been shown that at high 
energy collider a significant rate of events with double 
heavy quark pairs is expected. To compute all the fourth 
order ctg Feynman diagrams we have to consider the two 
different mechanisms leading to the same final state: the 
usual single parton scatterings and the double parton 
scatterings. At high energy in single parton interactions 
the partonic sub-process is dominated by the gluon-gluon 
fusion gg cccc while the double parton interactions are 
dominated by the process {g+g)-\-{g+g) {c+c) + {c+c) 
where the two distinct interactions occur in the same 
hadronic event. 

We give an estimate of the production cross section 
at high energy in the region of small transverse mo- 
menta where the multiple parton interactions provide the 



leading contribution to the cross section [9], [10]. We 
compute the production cross section of two c-quarks, 
ci,C2, very close in momentum space \pij — P2j\ < A, 
j = x,y,z, as a function of A. We consider the heavy 
quark production in the kinematical range of the LHCb 
= 14 TeV, 1.8 < 77 < 4.9), and for completeness for 
the ALICE (^i = 14 TeV, \r]\ < 0.9), Tevatron (^i = 
1.8 TeV, \y\ < 1) and RHIC = 200 GeV, \r]\ < 1.6) 
experiments; in the last case we calculate also the pro- 
duction cross sections in proton- nucleus interactions [? ] . 
The results are shown in Fig 1 One can notice that the 
cross section da/cPp at small A is almost rmiform and 
then it is approximately proportional to the momentum 
volume A^. 
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FIG. 1; Production cross section of two c quarks in momen- 
tum space A at LHC (LHCb and ALICE), at Tevatron and 
at RHIC. 

In the second step, the two c quarks join into a diquark. 
We assume simultaneous production of two independent 
c quarks with momenta Pi,P2- Since they appear wher- 
ever within the nucleon volume, we modulate their wave 
functions with a Gaus s ian prof ile with the "oscillator pa- 
rameter" B = y^73\/<~^^^~> = 0.69 fm corresponding 
to the nucleon rms radius 

where the normalisation factor A/"/? = 7r~^/^/3~^/^. Here 

fa is the average distance between two nucleons in the 
target nucleus for proton-nucleus experiment at RHIC 
and we use the value fa = fm or zero otherwise. 

We make an impulse approximation that this two- 
quark state is instantaneously transformed in any of the 



3 



eigcnstatcs of the two-quark Haniiltonian. Then the am- 
pUtude of the diquark formation M is equal to the over- 
lap between the two free quarks and the diquark with the 
same centre-of-mass motion. By approximating the di- 
quark wave function with a Gaussian with the oscillator 
parameter /? = 0.41 fm we get 




For the production cross section we take into account 
that da/d?p is practically constant and can be taken out 
of the integral 

where factors in front of the integral are due to the projec- 
tion on the colour and spin triplet states. If we insert the 
values of da/dp^ obtained from the Fig 1, we get cr « 27 
nb and 58 nb for LHCb and ALICE at LHC, a «:! 21 nb 
at Tevatron and cr « 4 nb and 63 nb at RHIC for proton- 
proton and proton-nucleus interaction, respectively. 

The last step of the Tec production is dressing of the 
heavy diquark. It either acquires one light quark to be- 
come the doubly-heavy baryon ecu, ccd or ccs, or two 
light antiquarks to become a tetraquark. With this we 
neglect the possible dissociation of the heavy diquark into 
a DD pair so the results are the upper estimate for the 
real T^c production. Assuming that the probability for 
dressing the cc diquark into the ccud tetraquark is 0.1 
[12], as pointed out at the beginning of the section, yields 
the production rate of the dimeson 20900, 9700, 600 and 1 
events/hour for LHC at luminosity 10'''^ cm~^s~'^. Teva- 
tron at luminosity 8 • lO'^^ cm^^s^^ and RHIC at Au 
luminosity 0.2 • 10^® cm~^s~^, respectively. 

III. STRUCTURE OF Tec 

The structure of the Tec tetraquark has been studied 
in rcf . [26] . We summarize here those features which are 
particularly relevant for the detection. 

There are two extreme spatial configurations of quarks 
in a tetraquark. The first configuration which we call 
atomic is similar to Ac, with a compact cc diquark instead 
of c, around which the two light antiquarks arc moving in 
a similar manner as in the Ag baryon. The second config- 
uration which we call molecular resembles deuteron, the 
two heavy quarks are well separated and the two light an- 
tiquarks are bound to them as if we had two almost free 
mesons. The atomic configuration is more likely to ap- 
pear in strongly bound tetraquarks while the molecular 
configuration can be expected in weakly bound systems. 

We present results using two different one-gluon ex- 
change potentials. The Bhaduri potential [18] quite suc- 
cessfully describes the spectroscopy of the meson, as well 



as baryon ground states. This is an important condi- 
tion since in the tetraquarks we have both quark-quark 
and quark-antiquark interactions. The ALl potential 
[19] slightly improves the meson spectra by introducing 
a mass-dependent smearing of the colour-magnetic term. 

We expand tetraquark wave function with Gaussinas 
of three sets of Jaccobi coordinates. In this basis we were 
able to reconstruct the wave functions of deeply bound 
tetraquarks as well as of two free mesons - the threshold 
state. This is important if one is searching for weakly 
bound tetraquarks with molecular structure. We found 
that the Tec is weakly bound for both the Bhaduri and 
ALl potential in contrast to the results of calculations in 
harmonic oscillator basis [21] where asymptotic channel 
cannot be accommodated as shown in Table I. 

TABLE I: Column 1: type of potential. Column 2: lowest 

moson-mcson threshold for a given potential in MeV, Column 
3: our results in MeV, Column 4: results in MeV of ref.[21], 
Column 5: mean distance in fm between two heavy quarks 

threshold our calc. Ref.[21] (fee) 
Bhaduri 3905.3 3904.7 3931 2.4 
ALl 3878.6 3875.9 1.6 



In Fig. 2 we present the probability densities pij for 
finding (anti)quark i and (anti)quark j at the interquark 
distance rij and the ratio of the projections on colour 
sextet state 1612634)0 and colour triplet state I3i2334)c 
where e.g. 

Ptr-\r) = (V.|3i2334)c(3i2334|c'5(r - 

Here particles 1 and 2 are the two heavy quarks c and 
particles 3 and 4 the light antiquarks u and d. The 
wave function between heavy quarks is broad and has 
an exponential tail ~ exp(— Kr) at large distances where 
K = ^y\Eh\Mred/^c, Eh is the binding energy of the 
system and Mred the reduced mass of the D and D* 
mesons. At small distances the dominant colour config- 
uration is 3i2334. Here we have a diquark-antidiquark 
structure and this region present about a third of the 
total probability while for r > 1 fm sextet colour con- 
figuration become larger. The ratio of these two con- 
figurations stabilise at 2, since here we have a molecu- 
lar structure of the two colour singlet mes ons which has 
in diquark antidiquark basis II13I24) = \/l/3l3i2334) -|- 
^2/3)612634) colour decomposition, while octet configu- 
ration I813824) = -v/273|3i2334) + \/T73|6i2634) is neg- 
ligible. 

Now we show that additional weak three-body inter- 
action can transform the molecular structure of the Tec 
tetraquark into atomic. For the radial part we take the 
simplest possible radial dependence - the smeared delta 
function of the coordinates of the three interacting par- 
ticles [20]. The colour factor in the two-body Bhaduri 
or ALl potential is proportional to the first (quadratic) 
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experimental values, this is actually a desirable feature. 



FIG. 2: Results for the ALl potential. Probability density of 
the two heavy quarks pec, of the two light antiquarks pqq and 
of a light antiquark and a heavy quaxk pqc in Tec as a function 
of the interquark distance. The ratio of the projection on 
colour sextet and colour triplet configurations is also shown. 
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FIG. 3: Results for the Bhaduti potential. Probability density 

between two c quarks pec in the Tec tctraquark as a function 
of interquark distance for three different values of the strength 
of the three-body potential. 



Casimir operator C^^^; C^^^ = A • A. It is then natural 
that we introduce in the three-body potential the second 
C(2) (cubic) Casimir operator C^^) = (P^'^Xa ■ \b ■ K> A 
deeper discussion of the properties that the colour de- 
pendent three-body interaction must fulfil can be found 
in [22-24]. 

In the baryon sector the three-body interaction was 
used to better reproduce the baryon ground state spec- 
troscopy [19]. A colour structure is there irrelevant since 
there is only one colour singlet state and thus the colour 
factor is just a constant which can be included into the 
strength of the potential. In tetraquarks the situation is 
different since there are two colour singlet states: 812334 
and 612634 (or I13I24 and 813824 after recoupling). The 
three-body force operates differently on these two states 
and one can anticipate that in the case of the weak bind- 
ing it can produce large changes in the structure of the 
tetraquark. This cannot be otherwise produced simply 
by reparameterization of the two-body potential, so the 
weakly bound tetraquarks are a very important labora- 
tory for studying the effect of such an interaction. 

A drastic change in the width of the probability den- 
sity can already be seen for strength Uq = -20 MeV (Fig. 
3), where the binding energy of the tetraquark becomes - 
14 MeV for Bhaduri potential. Here Tec loses the molecu- 
lar structure, the triplet-triplet colour configurations be- 
come dominant and the Tec tetraquark becomes similar 
to Thb- In the baryon sector such an interaction would 
merely lower the states by about Uq so it would have no 
dramatic effect nor would it spoil the fit to experimen- 
tal data. Since the predicted energies of ground state 
baryons for the Bhaduri and ALl potential are above the 



IV. DETECTION 

In order to identify a weakly bound Tec tetraquark we 
have to distinguish the pion or photon emitted by the D* 
meson bound inside the tetraquark from the one resulting 
from free D* meson decay. We can exploit the fact that 
the phase space for £)* D + n decay is very small. This 
has a strong impact on the branching ratio between radia- 
tive and hadronic decay. Since the D* meson inside the 
tetraquark with molecular structure is not significantly 
influenced by the other D meson in the tetraquark, we 
expect that the partial width for the magnetic dipole Ml 
transition would be very close to the width of the free me- 
son while the width for hadronic D* D + Tr decay will 
decrease with stronger binding and will become energeti- 
cally forbidden below the D + w threshold. The hadronic 
decay of the Tec tetraquark is a three-body decay which 
is commonly represented by the Dalitz plot. 

If the Tec tetraquark is below the D + D* threshold 
but above the D + D + j and D + D + ir, as was the case 
in our nonrelativistic potential models, the partial decay 
rate for the Tec ^D-|-D-|-7r is given by 
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where particles 1 and 2 are two final D mesons and par- 
ticle 3 is a TT emerging from the decaying tetraquark. 
Here mfj = {pd + Pd)"^ and 77133 = (pd + PttY and M 
is the mass of the tetraquark. Since the total masses of 
the D* + D and 2D + tt are so close there is a strong 
isospin violation in the decay which cannot be repro- 
duced with the Bhaduri or ALl potential where the D* 
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FIG. 4: Dalitz plot for four different values of the binding 
energy ^b- 
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FIG. 5: Width of the Tec tctraquark obtained with integration 
of the Dalitz plot J dmi2dm23/M^ where M = mo + mo* + 
El, is the mass of the tetraquark. The width at S;, = is 
normalized to be equal to the width of the free D* meson 
decay. 

and the D isospin doublets are degenerate. We shall not 

try to modify the interaction to accommodate the de- 
pendence of the decay on the isospin of the particles, 
but we shall rather work with the experimental masses 
taken from the PDG [25] where we see that mD*+ — 
m-£,+ — m„o = 5.6 ± 0.1 MeV, m,-£,,a — m,-£,o — m^a = 
7.1±0.1MeV, mD.+ - moo -m^+ = 5.87± 0.02 MeV. 
The allowed region of integration over dmf 2 and dm^s for 
three di fferen t binding energies is plotted in Fig. 4. If we 
assume is constant, which is very plausible in our 

case, the allowed region will be uniformly populated with 



experimental events so that the measured partial decay 
rate T will be proportional to the kinematically allowed 
area from Fig. 4. This is shown in Fig. 5, where we as- 
sumed that for the molecular state, the width of the T^c 
tetraquark with zero binding energy would be the same 
as the width of the free D* meson. 

Let us now consider also the posibility that Tec is not a 
bound DD* state but a resonant state above the D + D* 
threshold. Then if the resonance is situated near the 
threshold, there will be a significant fraction of hadronic 
Tcc^ D + D + w decays beside the Tcc^ D + D* decay. 
This region of positive binding energy is also presented 
in Fig. 4 and Fig. 5. 

In order to estimate the decay width we make a 
comparison with charmonium. The charmonium state 
^-(3770) has the width of 25.3 ± 2.9 MeV and is 36 MeV 
above DD threshold, which is also the dominant decay 
mode. Let us assume, that the Tec tetraquark resonant 
state, which would be 36 MeV above the D + D* thresh- 
old would have the same partial width for the decay into 
D and D* meson The area of the integrated Dalitz plot 
for the binding energy i?b = +36 MeV is then 37 times 
larger then at the threshold Ey, = 0. Since the experi- 
mental width for £>* ^ Dtt is 96 ± 4 ± 22 keV [25] we 
can estimate that the decay width for the Tec — * DnD 
three-body decay would be F ~ 37 • 96 keV = 3.6 MeV. 
So we expect about 15% direct Tec D + D + n decays. 
The Dalitz plot would not be uniformly populated but 
there will be a strong band where 77123 = ™d* reflecting 
the appearance of the Tec —* DD* — » DnD decay chain. 
In this estimation we have neglected the interference be- 
tween these two decay mechanism, since the width of the 
D* is three orders of magnitude smaller then the width 
of the tetraquark. 



V. CONCLUSION 

We have shown that the Tec tetraquark production is 
comparable to double charm baryon production (possibly 
10%). Therefore they may be seen in SELEX if statis- 
tics is improved. Similarly, it is comparable to prompt 
Jf-ipcc production which is reasonably abundant in B- 
factories. In high energy colliders we may expect an op- 
timistic number of events due to double cc production 
via double two-gluon fusion (see Sect. 2). Therefore time 
has come to start the hunt! 

Regarding the detection of the Tee = DD* tetraquark 
we propose a nice opportunity the very small phase 
space of the D* ^ Dtt decay which is very sensitive to 
the binding energy of D* to D. One possibility would be 
to measure the branching ratio between the pionic and 
gamma decay of £>*. 
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